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(54) Determining pore pressure while drilling a borehole 

(57) Pore pr^sure Ppore of a subsurface formation is detemiined, while drilling a borehole, by measuring 
formation strength FS, resistivity RES and natural gamma ray activity GR. 

The measurement FS includes contributions of clay volume V^, non-clay mineral strength, and a 
porosity, the measurement RES includes contributions from a porosity, non-clay mineral volume Vgnt and pore 
fluid conductivity and the measurement GR includes contributions from clay and non-clay mineral, and the 
sum of Vci, Vgjft and Ogff is 1. 

The contribution to the measurement of clay volume Vd, non-clay mineral volume V^wi, clay porosity 
4>ciay and effective porosity 0eff Bre derived for the formation being drilled, the clay porosity Oday and effective 
porosity 4>eff being related by » ^clayt'l " ^sitt)- 

The pore pressure Ppore Is derived from the relationship 
*eff = *clayW)n-Vsiit{Z)]expI-aPpore<Z)L 

wherein 



b» 
Pb = 

g = 

*eff = 
*cl8y = 
Vsllt = 



and 



compaction coefficient at depth of interest 

density of shale 

density of pore fluid 

gravitational acceleration, 9.81 ms*^ 

effective porosity 

clay porosity at zero depth of burial 

volume of non-clay mineral (silt) at depth of Interest 



W 

N) 

N3 
CJI 

CJI 
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PORE PRESSURE EVALUATION 



The present invention relates to a method of formation evaluation for use in 
drilling hydrocarbon wells or the like in underground formadons. In particular, the 
method provides a means for evaluating formadon pore pressure from volumetric 
properties such as porosity in shale formations. 

In the rotary drilling of hydrocarbon wells, drilling fluid or "mud" is circulated 
into the well through the drill string and returns to the surface via the annulus between 
the drill string and the borehole wall, bringing any rock cuttings with it The mud also 
serves to exert a hydrostatic pressure on the formation to provide physical stability to 
the borehole and to balance any pressurised fluids in the formation being drilled so as 
to prevent their influx into the borehole. This latter function is important for weU 
safety since if the formation fluid is a gas, the effect of an influx mto the mud will be 
to reduce its density and hence the hydrostatic pressure exerted on the formation. This 
in turn will allow more fluid to flow into the well further reducing the density of the 
mud and if remedial action is not taken sufflcientiy quickly a blowout can occur. 
However, it is also undesirable to exert too high a hydrostatic pressure on the 
formation since this can cause fracturing and formation damage due to fluid loss into 
the formation which is detrimental to the ability to produce hydrocarbons such as oiL 

As will be appreciated from above, it is generally desirable to maintain the mud 

density such that the hydrostatic pressure balances the pressure of the fluids in the 

formation. Consequently, it is desirable to be able to evaluate this pressure, known as 

''pore pressure" while drilling is being conducted such that the density of the mud can 

be controlled continuously in an appropriate manner. Previously proposed methods of 

evaluating pore pressure, such as those described in US 4,833,914 and US 4,949,575 

(incorporated herein by reference) rely on measurements made during drilling using 

MWD technology and formation properties from formation evaluation while drilling 

or wireline logs. The measurements are interpreted to obtain a volumetric evaluation 

of the formation in the manner known as "GLOBAL" and described in US 4,338,664 

which includes the porosity of the formation of interest since this is required to derive 

the pore pressure. The basic relationship used in this prior method to derive the 
effective porosity of the formation is: 

*<ir = *«4^-^d^r (1) 

where 4>„„ is the porosity of sand with no clay present and is the volume 
fraction of wet clay in the formation. 
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This interpretation relies on a model of a shaley sand to interpret the 
measurements. It is an object of the present invention to provide a method based on 
an alternative model so as to allow evaluation in cases where the prior technique is 
less appropriate. In the present invention, the model underlying the method is one of a 
shale having a proportion of silt (non-clay) material and the effective porosity is given 
by the relationship: 



(2) 



where O^^^ is tiie effective porosity of 100% clay at the depth or effective pressure of 

interest and is the volume fraction of non-clay materials. 

In its broadest aspect, the present invention provides a metiiod for determining 
the volumetric composition of a shale from its measured resistivity, wherein the 
measured resistivity RES is integrated according to the relationship: 



wherein is the shale resistivity, is shale conductivity, is the pore 

fluid conductivity, m is the power law exponent, is the volume fraction of non- 
clay mineral, is the effective porosity of 100% clay at the effective pressure of 

interest and is related to the effective porosity of the shale by the relationship: 

In accordance witii an embodiment of the present invention, there is provided a 
method of determiiung the pore pressure of a subsurface formation while 

drilling a borehole, comprising measuring the formation strength FS, the resistivity 

RES, and the natural gamma ray radioactivity GR so as to derive the contribution to 
the measur^ents of clay volume V^, non-clay mineral volume V^, clay porosity 
^ciay effective porosity for the formation being drilled, wherein the clay 
porosity and the effective porosity are related according to the relationship 
= ^ciayi}- " ^«ar) wherein the measurement FS includes contributions of clay 
volume V^, non-clay mineral strength, and a porosity, the measurement RES includes 
contributions from a porosity, non-clay mineral volume and pore fluid 

conductivity and the measurement GR includes contributions from clay and non-clay 
mineral, and the sum of V^, and being 1, the pore pressure being 

derived from the relationship 
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wherein a = 7 — ^ — r-and 

& = compaction coefficient at depth of interest 
= density of shale 
= density of pore fluid 

g = gravitational acceleration, 9.81 ms'^ 

<E>^ = effective porosity 

*c^y = clay porosiQr at zero depth of burial 

=s volume of non-clay mineral (silt) at depth of interest 

The method according to die present invention is applicable to cases where the 
pore pressure or overpressure is related to die porosity of die formation being drilled. 

Porosity is typically obtained from a nuclear log interpretation or optionally a 
resistivity log inteq)retation. Other suitable measurements typically include neutron 

porosity and total density. 

In the present method, each measurement made of the formation is related to die 
constitution of the formation by an equation, often known as a tool response equation, 
which expresses the measurement obtained in terms of the volumetric components of 
die formation. The tool response equations, in combination with an equation which 
states that the sum of the volumes of all of the components is equal to one, are solved 
simultaneously by an incoherence minimisation technique to produce a volumetric 
analysis of the formation. The measured parameters are combined with other 
information which is known for die formation of interest eg compaction coefficient, 
typically obtained from tabulated data so as to enable the pore pressure to be 
determined. 

The revised interpretation of shale resistivity provided by the present invention 

is related to the volumetric composition of the shale. The equation is based (m a 
model of shale as a matrix of compacted clay, whose porosity is O^, containing 

isolated spherical particles of silt: 

y—=<^M. = ^fi^id^doTi^'^iaiV (3) 

where p,^ is shale resistivity, a,^ is shale conductivity, Gjj^ is die pore fluid 

conductivi^ and m is the power-law exponent 

The conductivity of the pore fluid is determined by first calculating its molarity 
which is then interpreted using tabulated ionic conductivity data and corrected for 
downhole temp^tuve. 



The simple method for calculation of total molarity is just to add the exchange 
ion molarity, calculated from CEC (cation exchange capacity for clay minerals) and 

porosity, to the deposition/equilibrium fluid molarity, and calculate conductivity on 
this basis. Using this route, the pore fluid molarity is given by: 

where Af^ is the depositional or equilibrium fluid molarity, CEQ and are 

respectively the CEC (equivalents/unit mass) and mass fraction of the i'th mineral in 
the clay matrix, and the mean density of the solid material in the clay matrix. In 
practice, is likely to be a free parameter which must be determined from 

appropriate cross plots. Its effect is, in many cases, small compared to that of the CEC 
term. 

A more complex but more conect route is to use the Donnan theory of 

equilibrium of membranes to calculate pore fluid molarity. In relatively fresh fluids, 

the simpler route is adequate. Either way, the conductivity of the pore fluid will be 

dominated by the motion of cations (since the clay contains a large number of fixed 

macro-anions) and the appropriate tabulated ionic conductivities which can be 

obtained from tests on representative samples are preferred rath^ than simple molar 

conductivities. If this data is limited, assumptions must be made about the partition of 

conductivity between cation and anion to obtain an estimated value. 

The tabulated values are quoted for a ten^rature of, say, Tq, and the value for a 

downhole temperature T can be approximated using an activation energy term ie, 



(5) 



The activation energy term E^ can be obtained from observing the variation in 

conductivity of fluid at differing temperatures and is known not to vary to any great 

degree over the scale of interest in this application. Consequendy, a single 
representative value of E^ can be obtained by e^qieriment and applied generally in a 

region of investigation. The downhole temperature T is often difficult to measure 
directiy but can be inferred from a determination of the geothermal gradient obtained 
from wells in the same region. Since the temperature T can have a significant effect, 
the best estimation of T is desirable. 

The only remaining unknown in the resistivity equation is the power-law 
exponent m. This is known to vary widi day mineralogy, but can be estimated over 
an interval from a cross-plot of 



log 



(6) 



Any other depth-varying data, log-derived or otherwise, can be incozporated 
into the *GLOBAL* processing method. For example, in addition to gamma-ray, 
formation strength and resistivity required for the method of the present invention and 
the pore pressure determination described in U.S. patent 4,949,575 the response 



where GR is total gamma-ray response, is the volume fraction of mineral i in the 
rock, and GR^ is the gamma-ray response of that mineral in the absence of the others 
(if spectral ganmia-ray data is available, then a similar set of mineral based response 
equations can be written for each energy window recorded by the tool, e.g. K, U and 
Th decay, which give a more detailed set of inputs to the optimization procedure); 



where FS is formation strength (derived from MWD drilling parameters and 
normalized for variations in mud and pore pressures), FS^ is the formation strength 
of the non-clay mineral, is the volume of wet clay in the formation, is the 

clay porosity, and the quantities with 'zero* subscripts are the values of these 
parameters, picked from cross-plots, at which the formation strength falls to zero; 



where NPHI is the total neutron porosis response, and PN^ and are respectively 
the neutron pcnosity response and the volume fraction in the formation of the Tth 
component of the formation (these components are pore fluid, the various silts and 
clays making iip the shale, and organic matt^); 



equations for neutron porosity and gamma density can also be used and can be 
applied usefully to the method of the present invention. The basic response equations 
are respectively: 



GR^^VJGR, 




RHOB^^ViRHOB^ 
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where RHOB is the total density response, and RHOB^ is the density response of the 
i'th component as above. 

The GLOBAL method allows diese and other available logs or similar data to 
be combined in the volumetrics determination, provided a suitable response equation 
exists. Logs which provide mineralogical discrimination such as spectral gamma-ray 
and photoelectric cross-section are particularly useful for the present method; for 
example, resolution of total clay into distinct clay minerals can give a better value for 
CEC of the rock, and hence make the resistivity equation, which is sensitive to this 
parameter, a more effective input in the volumetrics process. By placing a volumetric 
constraint on the formation responses, violation of the constraint causes the 
incoherence determined by GLOBAL to increase which in turn causes the 
minimisation to bring the volumetrical values back into the constrained levels. A 
continuity constraint can also be included to prevent wild fluctuations in calculated 
values in adjacent depths which can further improve the calculated results. 



Qaims 



A method of determining the pore pressure P^,^ of a subsurface formation 

while drilling a borehole, comprising deriving the contribution to the 
measurements of clay volume V^,, non-clay mineral volume V^, clay porosity 
O^i^ and effective porosity for the formation being drilled, wherein the 
clay porosity *^ and tfie effective porosity are related according to the 

relationship 

by measuring the formation strength FS, the resistivity RES, and the natural 
gamma ray radioactivity GR and wherein the measurement FS includes 
contributions of clay volume , non-clay mineral strength, and a porosity, the 

measurement RES includes contributions from a porosity, non-clay mineral 
volume V^fjf and pore fluid conductivity and the measurement GR includes 
contributions from clay and non-clay mineral, and the sum of V^, and 
being 1, the pore pressure being derived from the relationship 

=«»*»(0)[l-V^(Z)]exp[-aP^(Z)]^ 

wherein a = 7 — - — r" ^d 

b = compaction coefficient at depth of interest 

s density of shale 
p^= density of pore fluid 

g = gravitational acceleration, 9.81 ms*^ 
s= effective porosi^^ 

^etoy = cl^y porosity at zero depth of burial 

= volume of non-clay mineral (silt) at depth of interest 

A method as claimed in claim 1, wherein die formation strength {FS) is derived 
from measurements of downhole bit torque (TOR) and downhole weight on bit 
(WOB) and wh^ said formation strength is corrected for the effects of bit 
wear. 



A method as claimed in claim 1 or 2, wherein 

1--Jld SSL. 



FS = FS, 



where FS is formation strength, FS^ is the formation strength of the non-clay 
mineral, is the volume of wet clay in the formation, is the clay 

porosity, and the quantities with 'zero* subscripts are the values of the 
parameters at which the formation strength falls to zero. 

4 A method as claimed in any preceding claim, wherein 

GR = ^VfiR, 

where GR is total gamma-ray response, Vj is the volume fraction of mineral i in 
the rock, and GR^ is the gamma-ray response of that mineral in the absence of 

any others. 

5 A method as claimed in any preceding claim, wherein 

PshaU 

where is shale resistivity, (T^ is shale conductivity, cr^ is the pore 
fluid conductivity and m is the power-law exponent 

6 A method as claimed in claim S, wherein, for a temperature (T)* 



where (Tjtdd{TQ) is the pore fluid conductivity at a temperature T^, is an 
activation energy temi mdR is the gas constant 

7 A method as claimed in claim S or 6, wherein the power law exponent n is 
obtained from a cross-plot of 



8 A method as claimed in any preceding claim further including measurement of 
one or more of neutron porosity, gamma ray density, sonic travel time deep 
induction resistivity and resistivity at bit 

9 A method as claimed in any preceding claim, wherein FS is derived as a 
function of rate of bit rotation, bit efficiency, gouging component of bit torque, 
rate of penetration, bit diameter and drilling mud weight 
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A method of determining the volumetric composition of a shale comprising 
relating the effective prorosity of the shale to the porosity of 100% clay at 

the pressure of interest O^,^ according to die relationship: 

wherein is the volume fraction of non-clay mineral and measuring the 
resistivity RES of the shale which is interpreted according to the relationship: 

wherein is die shale resistivity, a,^ is shale conductivity, (T^ is the 
pore fluid conductivity and m is the power law exponent 
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